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Mean Field Dynamics in Non-Abelian Plasmas from Classical Transport Theory 
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Based on classical transport theory, we present a general set of covariant equations describing the 
dynamics of mean fields and their statistical fluctuations in a non-Abelian plasma in or out-of- 
equilibrium. A procedure to obtain the collision integrals for the Boltzmann equation from the 
microscopic theory is described. As an application, we study a hot non-Abelian plasma close to 
equilibrium, where the fluctuations are integrated out explicitly. For soft fields, and at logarithmic 
accuracy, we obtain Bodeker's effective theory. 
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The dynamics of mean fields in non-Abelian plasmas 
is at the basis for an understanding of many proper- 
ties of the early universe. This concerns the quark-gluon 
plasma which is believed to be formed at high temper- 
ature and/or large chemical potential, or the plasma 
formed in the symmetric phase of the electroweak the- 
ory. The energy densities required to form a quark-gluon 
plasma might be reached in the coming generation of 
heavy ion colliders. A theoretical framework to describe 
the dynamics in non-Abelian plasmas in or out of equi- 
librium is thus mandatory (see |l| for an early review). 

The dynamics of soft fields with momenta about g^T 
is dominated by non-perturbative phenomena, even for 
small gauge coupling g <^ \ and close to equilibrium 
H]. Here, some progress has been achieved recently 
by Bodeker ||], who derived an effective theory for the 
soft modes after integrating-out the perturbative physi- 
cal scales T and gT from the field theory. 

In this Letter, we study mean field dynamics in the 
context of classical transport theory ji) . In the close-to- 
equilibrium plasma, classical transport theory is known 
to describe correctly the physics around the scale gT |y , 
in full agreement with the results obtained from quan- 
tum field theory g]. We will show how also the physics 
at scales g^T, and at logarithmic accuracy, can be under- 
stood in the same kinetic language. Our ultimate goal 
is providing a well-defined prescription to treat a num- 
ber of interesting problems including out-of-equilibrium 
phenomena. Here, we present the main results of our 
analysis. This includes the derivation of a general set of 
dynamical equations for the mean fields and their fluctu- 
ations, and a prescription for obtaining collision integrals 
for Boltzmann equations. An explicit application for hot 
non-Abelian plasmas close to equilibrium is given as well. 
We leave for a longer publication a more extensive dis- 
cussion and further details on the computations M] . 

We begin with the derivation of a set of equations 
that describe both the dynamics of fluctuations and of 
the mean fields in a non-Abelian plasma. We follow 
the same philosophy and methods which have been used 
longly in kinetic theory to study Abelian plasmas (see 
g]). Our equations can be seen as the generalization to 
non-Abelian plasmas. 



The starting point is to consider gluons (or electrons, 
quarks) as charged point particles moving on a world line 
according to their classical equations of motion. They in- 
teract via a classical chromo-electromagnetic field. We 
consider a system of these particles carrying a non- 
Abelian color charge Q°' , where the color index runs from 
a = 1 to N'^ — 1 for a SU(7V) gauge group. 

Within a microscopic description, the trajectories in 
phase space, and therefore the one-particle phase space 
density f{x,p,Q), are known exactly. The classical tra- 
jectories x(t),p{t) and Q{t) are solutions of the Wong 
equations 
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Here, A" denotes the microscopic vector gauge field, 
F;,[A] = d^^Al - d^Al + gfabcA'^^Al the corresponding 
microscopic field strength, and /"''^ are the structure con- 
stants of SU(iV). We set c = ?i = /cb = 1, and work in 
natural units. Effects of spin will be neglected. 

The microscopic phase space density obeys Liouville's 
theorem df/dr = H. We can write it using (nl) as M 



[d^ - gr'^A^Q^dQ ~ gQaF;:,d;] f = 



(2a) 



with dP = d/dpf" 



and d^ = 9/9(3°. This equation is 
gauge invariant, with / transforming as a scalar B, and 
is completed with the Yang-Mills equation. 



D^pt^'^ix) = r{x) 



(2b) 



In the adjoint, the covariant derivative is given by 
D'^''[A] = d^S"" + gf^^A^. In a self-consistent picture 
the current Jl^{x) is due to the particles themselves, thus 



Ji:ix)=9 Y. fdPdQp''Qaf{x,p,Q). (2c) 



(Throughout, we shall omit a species or helicity index 
on /, and the sum over species and helicities will, in 



the sequel, not be given explicitly.) Physical constraints 
are enforced through the phase space volume element 
dP = d^p 2Q(po)5{p^ - m2)/(27r)3, while dQ contains ^- 
functions imposing the group Casimirs (see ||5|| for their 
definition). The covariant conservation of the current 
( Pq) is shown using (|2^) g . 

If the system under study contains a large number of 
particles it is impossible to follow their individual trajec- 
tories in phase space. Thus, / can no longer be consid- 
ered a deterministic quantity and one has to switch to a 
statistical description, taking statistical averages (. . .} of 
all microscopic quantities. We write 

Al ^A^+al, f = f + Sf , J„^ = J^ + SJ^ , (3) 

where the quantities with a bar denote the mean values, 
e.g. f = (/) and A = {A), while the mean values of 
fluctuations vanish, {Sf) = and (a) ~ 0. We also split 



pa ^ pa I ra 






(4a) 
(4b) 



with D = D[A] and F_ = F[A]. Note, that the mean 
field strength (i^^^) = F^^ + g f"^" {a^al) due to the non- 
Abelian nature of the theory. 

Let us take a statistical average of (g) to find the ki- 
netic equations for the mean values, 

P^ (D^ - gQaF^.d;) f = (r;) + (C) , (5a) 

D^F^-- + ( JIJ = r . (5b) 

In (1^ we used [9^ -^/''^QcA^^^]/ = D^f. The tunc- 
tions 77,^ and J^^^ are of second and higher order in the 
fiuctuations and read 



gQaP''d;f",,6f , 



^:^f-i' £abcr^c 
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r'^Q^idQa^Sf +galald;D 



(6a) 

(6b) 



jr ^ .9/''^ [^L«^M«c + Sada,,, fr] • (6c) 

The corresponding equations for the fluctuations are ob- 
tained by subtracting (ph from (0). The result is 

p^ [D^ - gQaF;:,d;] Sf - gp^at^^f-'^Qcd^f 
-gQa[D,,a, - D,a^]y^^dPf = ,7 -^ ^ - (?? + (7a) 

+Jn:'^ - (^fluo') = SJ'"'" ■ (7b) 
A number of comments are in order. 

1. The equations (|5|) and (|^) are exact, no approxima- 
tions have been made. In particular, they are also 
valid in out-of-equilibrium situations. 

2. The equations (g|) and (0) are consistent with gauge 
invariance. They are covariant under the mean gauge 
field symmetry SeA"^ = {D^eY and S^a^ = 5/"'"=a|^e'=, 
in analogy to the background field formalism [nO| . This 
establishes the compatibility of the statistical average 
with the gauge transformations of the mean field. We 
postpone a careful and detailed discussion to . 



3. The microscopic current conservation implies 

{D^.P)^+gfabc{alS.r-^)=0. (8) 



This is automatically consistent with (pH) , provided J 
and SJ are solutions of (|^) and (|^) . (A similar equation 
holds for the fluctuation flelds.) Note that the validity 
of (pf) turns into a non-trivial consistency check for 
approximate solutions. 

4. The functions {rj) and (^) can be considered as the 
effective collision integrals of the Boltzmann equation 
(pa|). In our formalism the collision integrals arise as 
correlators of statistical fluctuations. The fluctuations 
of the gauge flelds cause random changes in the motion 
of particles, and thus, they can be viewed as having the 
same effects as collisions. This can be seen as a deriva- 
tion of collision integrals from the microscopic theory. 
Note also that the current induced by the fluctuations 
of the gauge field (Jtiuc) is a purely non-Abelian effect. 

5. A general procedure for integrating-out the fluctua- 
tions amounts to flrst solve their dynamics (|7|) in the 
background of mean fields. In general, this is a difficult 
task, in particular due to the non- linear terms in (M). 
The (explicit) solution is then inserted into (g). The 
back-coupling of the fluctuations to the mean fields 
is finally obtained after taking the statistical average 
of the functions (|g), and yields the effective collision 
integrals and the induced current in (0). 

6. The set of equations ^ and (H) reproduces the known 
set of kinetic equations for Abelian plasmas in the cor- 
responding limit PI, in which only the collision in- 
tegral (77) survives. The Abelian counterpart of (77) 
can be expressed as the Balescu-Lcnard collision in- 
tegral |g[. One can then proof in a rigorous way the 
correspondence between fluctuations and collisions in 
the Abelian plasmas mentioned above. (An analogous 
derivation of collision integrals for Wigner functions 
can be found in flit], see also p2[.) 

7. Neglecting all fluctuations reduces (H) to the well- 
known (non-Abelian) Vlasov equations. 

This terminates the derivation and discussion of the 
basic set of equations. 

To put the method to work we will specialize our analy- 
sis to hot non-Abelian plasmas close to equilibrium, with 
the gauge coupling g <C 1. This allows us to perform sev- 
eral approximations. We will consider small fluctuations, 
neglecting in (^^ and ( pq ) the terms cubic in the fluc- 
tuations. This is interpreted as neglecting effective three 
body collisions versus binary ones. In the same spirit, 
we employ the second-moment approximation for the dy- 
namics of the fluctuations [p|, setting 1] = {rj}, S, — (^) 
and Jfiua = {Jiiuc) in (0). This linearizes the dynamics of 
the fluctuations and can be interpreted as neglecting the 
influence of collisions on the dynamics of the fluctuations. 



Finally, the term containing the mean field strength in 
( [7a| ) is negligible compared to the remaining terms and 
will be omitted, as long as glF^^l/niD ^ T, with mu 
the Debye mass Q. 

We study the mean distribution function f(x,p, Q) = 
I°'^'{Po) + 9l^'^\x,p,Q). In the strictly classical ap- 
proach, the relativistic Maxwell distribution at equlib- 
rium is used for any species of particles. Here, we 
consider only massless particles in the adjoint, with 
J dQ QaQb — N Sab- For particles in the fundamental 
one has / dQ QaQb = | Sab instead. Solving (||) for van- 
ishing fluctuations in this approximation gives the infinite 
set of non- Abelian hard thermal loops MM . 

We now include small statistical fluctuations 6f around 
/ and re- write the approximations to (|5|) and (|7|) in terms 
of current densities and their fluctuations. Consider the 
current densities 

J'a,-aSx,p)^gpP fdQQa,---QaJ{x,p,Q), (9a) 



(9b) 



JS,-aSx,v) ^ dPJP^...aJx,p) . 



The measure dP only integrates over the radial compo- 
nents, dP = dPdn/4n, and w^ = (l,v) with v^ = 1. 
The current (^) is obtained performing the remaining 
angle integration J{x) = J ^^{x, v). From now on we 
will omit the arguments of the current density J, unless 
necessary to avoid confusion. 

After multiplying (|5^) by gQaP^ Ipq, summing over the 
two helicities, and integrating over dPdQ, we obtain for 
the mean current density at leading order in g 

v^D^jP + mlvPv'^F^a = (r^P) + i^) , (10a) 

D^F>'- + {JU = J" , (10b) 



with the Debye mass mj^ 
and 
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C'a^-gfabcVPa'l.SJ-'P, (lib) 

Jp^^^^ = gf^'^lD^^'^a^aP + 5'"^al{DPaP - ^^a^)'} . (lie) 
For the fluctuations we find 



-glabcV^aU'^P , 



^^^ , (12a) 

; ^bd + fmbdSac) J cd^ (12b) 

{&a^ - DP{Da)\ ^ + 2gfabcFrac^, = SJP . (12c) 



[v''D^SJP]^^^gv'^a";{f„ 



We solve the equations for the fiuctuations (Oh with 
an initial boundary condition for Sf, and a^(i = 0) = 0. 
Exact solutions to (12a) and ( |12b| ) can be obtained m. 
The current fluctuation SJa reads, for xq = t > 0, 



SJPix, v) = Uab{x, Xt) 5jP{xt,v) 

dTUab{x,Xr)\gfbdcV^a'lXxT)Jc{xT,v) 

+ mlvPy'^ {D,,ao - 5oa^)' (x,)} . (13) 

We introduced Xr = x — vt and the parallel transporter 
Uab, obeying vPD^Uab{x,y)\y=x^ = 0. In order to solve 
( 12c ) for a^ , we make a double expansion in both A and 
J, using Uab = Sab + 0{gA). We denote by a'") the 
term containing a total of n powers in the mean fields A 
and/or J. 

For our purposes, it will be sufficient to consider the 
zeroth order term in A, and the zeroth and first order 
terms in J. Using the one-sided Fourier transform 1^, 
we find 



<n(fc) = 



1 



-fc2 + Ut 
T{1) n,\ -gfabc 



An 



dn^ S J,'^a{t = 0,k,v) 

—ik-v 

dn^ 1 

1 ^ 

An —ik-v 



^T(k)/ 



F + Ht 



(14a) 



(14b) 



in the gauge k • a — 0. The function IlT(k) is the 
transverse polarization tensor of the plasma, Pf^(k) — 
5ij — kikj/'k'^ the transverse projector, and a[ = P^aj- 
Retarded boundary conditions are assumed above, with 
the prescription kg —> kg + iO^ . 

With the above, we can express all fluctuations in 
terms of initial conditions and the mean flelds. Follow- 
ing [pi the statistical average over initial conditions can 
be deduced m and flnally expressed (for each species or 
helicity index) as 



5(3) (p _ p')5(Q _ Q') J^p) + i,^j,.p.QQ, 



(15) 



The function 5{Q — Q') is, apart from a (representation 
dependent) normalization constant, a iV^ — N dimen- 
sional (5- function over the proper set of Darboux variables 
related to the color charges 0. The second term in (nS) 
is the Fourier transform of a smooth function that van- 
ishes at large distances. The above statistical average is 
all we need to evaluate the collision integrals. 

For the remaining part we will concentrate on the dy- 
namics of mean fields with typical momenta around gmjj . 
When computing the related collision integrals, we will 
find logarithmic divergences, cut-off in the infrared by the 
inverse collision time. We employ the leading logarithmic 
approximation, assuming ln(l/g) ^ 1 while neglecting 
all sub-leading (though finite) terms. 

We find that the induced current (JfLc) vanishes, as do 
the fluctuation integrals (t]^^)} and (■f'-*'-'). The vanishing 
of (fy*-"-*) is consistent with the fact that in the Abelian 
limit the counterpart of (rj) vanishes at equilibrium H. 



In the same spirit we evaluate the colhsion integrals con- 
taining one J field. Consider 



/•oo 

+gfcdev'' / dTjp,e{xr,v) ^a^°|,(a;) a';^l{xr)jj, (16) 





which simplifies, at logarithmic accuracy, to 

^2 
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e;,l{x,v)) =-^NT\nil/g)x 
dilv' 



V, 
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I{v,v')JS{x,v'), 



TT ^1 - (V V')^ 



(17) 

(18) 



The above expression has been obtained first in 0, and 
reproduces the collision integral considered in the Boltz- 
mann equation of p3| . 

We verified that the leading logarithmic solution is con- 
sistent with gauge invariance. Evaluating th e cor relator 
in (0) yields D^J'^ = 0, in accordance with (10b) in the 
present approximation. 

Following Bodeker, one can now estimate ,P from ( |lOaD 
to obtain for (10b) 



{D.Fna^aEl + ul, o 



inmjy 



3Ng^T\ii{l/g) 



(19) 



This is the result of Q. The coefficient a represents the 
color conductivity and has been discussed in ilq,H . The 
white noise ly has its origin in the fluctuations of the 
transverse part of ^^°^ PMI- We obtain to leading order 



K{x),yliy))^2TaS^^6ab6^^H^-y) 



(20) 



in accordance with the fluctuation-dissipation theorem 
(FDT). Note also that the classical Debye mass differs 
from the quantum one. 

In order to go beyond classical transport theory we ex- 
pand about the bosonic (fermionic) quantum-statistical 
equilibrium distribution function /+ (/_). For gluons in 
the adjoint, the Debye mass obtains as m|, = g^NT^/3. 
The FDT is obeyed as well, if / in ( [l5| ) is replaced by 
/±(1 ± f±). (This should however be derived in a sim- 
ilar way as ( |l5| ) from the microscopic theory [Q.) Also, 
the quantum collision integrals are obtained with the cor- 
rect statistical factors m. It is interesting to note that 
all quantum modifications are contained in the implicit 
change of m]j. 

This terminates the explicit derivation, in the leading 
logarithmic approximation, of the collision integral and 
the dynamical equations for the soft fields from classical 
transport theory. 

Summarizing, we have given a prescription to derive 
mean gauge field equations from classical transport the- 
ory. This includes a recipe to obtain effective (classical or 



quantum) collision integrals from the microscopic theory. 
The approach is in accordance with gauge invariance. In 
a close-to-equilibrium plasma and for small gauge cou- 
pling, we reproduce Bodeker's effective theory. 

The last part of our analysis can straightforwardly be 
generalized in order to obtain explicit expressions for the 
collision integrals not only for the soft momentum re- 
gion. Another interesting open problem is using the same 
methods for out-of-equilibrium situations. Based on the 
evaluation of collision integrals for Abelian plasmas out 
of equilibrium B , we should find the Coulomb logarithm 
changing drastically the mean non-Abelian gauge field 
equations. 

It remains remarkable that classical transport theory is 
efficient enough as to describe not only the non-Abelian 
dynamics of semi-hard modes with momenta around mo , 
but as well the non-perturbative dynamics of soft glu- 
ons at leading logarithmic order. This establishes a link 
even beyond the one- loop level between our approach and 
a complete quantum field theoretical treatment, whose 
deeper structure is waiting for being uncovered |1^ . 
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